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ABSTRACT 

The problem of applying the theory of systems analysis 
to the marine propulsion plant Is Investigated In this thesis. 
It Is believed by the author that the use of this theory can 
be of considerable assistance to the naval engineer In giving 
him an tinders tending of the fundamental behavior of the system 
and will provide a scientific basis for design criteria and 
operational policies. 

There are five basic areas of Investigation covered In 
the thesis which all contribute to the organization of the 
overall method of accomplishing an operational analysis of a 
complex mechanical system. The following Is a concise digest 
of each of the sections. 

1. The derivations of the basic mathematical formulas 
for reliability determination are given. The concept of com- 
ponent failure distribution Is presented and the exponential, 
normal, and Welbull failure models are discussed. 

2. The Importance of defining the objectives of the 
analysis clearly and completely Is shown. The Interpretation 
of component operating conditions and the resolving of 
success and failure standards are explained and an outline 
for the formulation of the necessary Input data required to 
undertake the mathematical part of the analysis Is offered. 

A hypothetical combined diesel -gas turbine propulsion plant 
Is used to Illustrate the proper procedures for performing a 
reliability analysis. 

3. The mathematical procedures for determining 
system reliability by the correlation of the Individual 
system components behavior are explained. In addition to 
presenting methods for solving serial, parallel, and Bayesian 
systems, formulas are derived for the reliability analysis of 
periodically operating components and of standby units. 
Representative values are postulated for the dlesel-gas 
turbine plant and reliability curves are plotted which demon- 
strate the effect of component wear-out, reliability as a 
function of time, and the Individual effect of specific com- 
ponents on the overall system reliability. 
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4 , Two techniques for ascertaining maintenance 
schedules are formulated. The first technique yields the 
optimum time to perform preventive maintenance in order to 
maximize system availability or minimize expected repair 
costs. The second procedure deduces by dynamic programming 
an optimal decision policy on whether or not to repair the 
system dependent upon the nmber of trips made since the last 
repair and the number of trips still to be accomplished. 

5 . A brief investigation is made into the problem 
of components subject to the interaction of their failure 
distributions. Pour possible interaction fimctlons of 
Instant failure rate are proposed and analyzed as to their 
effect on reliability and mean time to failure. An example 
is used to show how important this interaction effect can be 
when it is time dependent, such as would be caused by gradual 
wear. 



In conclusion this thesis presents a procedure for ana- 
lyzing a complex mechanical system using hypothetical 
examples to better illustrate the utility of the results to 
the naval engineer. 



Thesis Supervisor: Ernst G. Frankel 

Titles Assistant Professor of Naval Architecture 
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NOMENCLATURE 


Symbol 


Meaning 


R(t) 


Reliability 


p(t) 


Instant failure rate 


f(t) 


Probability density function of 
failure 


P(t), Q(t) 


Probability distribution of 
failiire 


T^, MTBP 


Mean time before failure 


t, T 


Time 




Exponential instant failure rate 


C 


Constant 



Specific nomenclature for diesel -gas turbine reliability analysis 



<^1 


p(t) for diesel engine 




p(t) for gas turbine 


^3 


p(t) for heat exchange 


C4 


p(t) for reduction gear 


s 

^6 


p(t) for lube oil pump 

p(t) for standby lube oil pump 


^7 


p(t) for fresh water pump 


% 


p(t) for salt water pump 


N 


Mean time of non use of the gas 
turbine 


^2 


Probability that the gas turbine 



will not be in use at time t 



INTRODUCTION 



Operational Analysis and Reliability Engineering 

During and shortly after the Second World War, due to 
the growing complexity of all types of operations In 
general, a branch of mathematics termed operations analysis 
became prominent. Probability theory forms the basis for 
this field of science. If a definition is desired. It may 
be said the operations analysis is a scientific study of 
operations to better londerstand their behayior for the 
twofold purpose of predicting their future results due to 
changes in the system and controlling the operation to im- 
prove its result. This thesis is concerned with the use of 
operations analysis methods to solve problons in the field 
of reliability engineering. 

In any system, mechanical or- otherwise, one of the most 
Important parameters determining its overall v^ilue is the 
reliability that the system can be expected to have during 
its operation. Although the exact definition of reliability 
is usually dependent upon who is the definer, the accepted' 
definition endorsed by the American Society for Quality 
Control [ 1 ] states that reliability is the probability of a 
device operating within specified limits for the time and 
operating conditions imposed on the device. Operations 
analysis techniques are used in reliability engineering for 
the purpose of predicting reliability and optimizing the 
system from a reliability standpoint. While reliability 
analysis is a relatively new;, subject which has generally 
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only been applied to the electronics and missile field, it 
has the potential of being extremely useful in the vinder- 
standing and improving of all types of mechanical systems 
including marine propulsion plants. 

Importance of Reliability Analysis 

The ability to be able to predict reliability is tanta- 
mount to understanding the underlying causes of \mrelia- 
bllity. The consequences of unreliability are numerous j the 
most obvious ones being cost, lost time, and danger to 
personnel. The cost Involved not only is that of the device 
to be replaced or repaired, but also the cost of the mainte- 
nance in material, manpower, and training with the added 
expense of having to continuously monitor the equipment for 
possible failures. The value of the analysis is that in 
choosing some specific system to perform a task, a relia- 
bility analysis of the system will furnish data on which to 
base a decision. For Instance, the Armed Services mainly 
places the emphasis for selecting a component or system on 
the low bid price. Normally, reliability is only given 
token consideration in the decision. The question arises as 
to whether a reliability investigation of the device should 
be made since isn't it possible that it would be advisable 
to pay more initially for a reliable system that in the long 
run will give better availability and less maintenance costs 
than a cheap unreliable system? Conversely, the possibility 

i , 

arises as to whether in some applications too much relia- 
bility is required of a device. Certainly the adage that 
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"a chain is only as strong as its weakest link" applies in 
a system and, for example, designing a valve with a safety 
factor of 5 when the piping only has a safety factor of 2 
is a waste of time and money » It is therefore evident that 
reliability apportionment should enter into any preliminary 
system design. 

Another problem inherent to this day and age is that of 
complexity. The major problems to be solved have become 
more intricate as time passes and it is only natural that 
there is an increase in the complexity of the equipment 
designed to meet these problems. The more complex equipment 
is, the greater number of components are required, and the 
more likely it is that one of these components will fail. 
Thus, an increase in complexity leads to a decrease in 
reliability. Of what value is complex equipment with theo- 
retically high performance if failures of its components 
keep the piece under repair most of the time. It can 
therefore be seen that engineering improvements for relia- 
bility will in many cases improve operability, maintaina- 
bility and productivity with an overall net gain in system 
worth . 

Other System Analysis Results 

Results other than reliability may also be obtained from 
an operations analysis of a system. The question of main- 
tainability of the system should normally be investigated 
unless for some special reason the system cannot be main- 
tained during its operating cycle, i.e, missiles, torpedoes. 




r 
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etc. A study of maintainability will give answers to such 
questions asi 

1. Is preventive maintenance worthwhile for the 
particular system? 

2. If worthwhile at what time intervals should it be 
scheduled? 

3. What is the optimum efficiency that can be expected 
from a system preventively maintained? 

In conjunction with preventive maintenance scheduling, 
the renewal theory of operations analysis [ 2> ] will give the 
expected number of replacements for a component over its 
lifetime, t. This gives a good basis for making an estimate 
of the number of spare parts to be stocked. 

Some other results that are obtainable arej 

1, The effect of redundancy on reliability. 

2, The system availability that can be expected. 

3, Optimization procedures for reasons of reliability 
or cost. 

4, Confidence levels that can be applied to the 
results of the analysis, 

5, Information for the determination of equipment 
quality control specifications. 

It is seen that by the methods of operations analysis, 
a deep Insight into the behavior of a system can be obtained 
which will produce not only a better understanding of the 
system but also give a firm basis for making decisions to 
improve its overall value. 
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Present State of the Art 

At the present time the state of the art of operational 
analysis of mechanical systems is very sketchy. Many mathe- 
matical models have been devised of hypothetical systems 
which have been analyzed thoroughly but the analysis of real 
systems has only been accomplished in the electronics and 
missile fields and even there only in limited situations. 
There are two main reasons for this lack of progress. First 
and foremost is the fact that the value of systems analysis 
and especially reliability engineering has just in the past 
few years come into prominence. Now that the mathematical 
tools have been developed and the value realized, the majori- 
ty of large companies are beginning to organize reliability 
analysis departments. The second factor and possibly the 
hardest to solve is the cost and time involved in running 
the testing programs necessary for the analysis. In order 
to obtain any results from an analysis, fail-ure tests must 
be run for the devices under Investigation and reams of data 
must be evaluated. Although this statistical aspect of 
operations analysis will not be discussed in this thesis, 
it should be mentioned that there are sampling procedures 
that can be used effectively to make data correlation ef- 
ficient. Inherently, the cost of initiating a reliability 
program will be high but just as inherently is the fact that 
in the long run, the program will more than pay for itself 



in its results 



I 
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Thesls Ob.lectlve 

The object of this thesis is to show how the methods of 
operations analysis can be employed in the investigation of 
a marine propulsion plant and of what worth such an investi- 
gation will be. The overall study of this question will be 
divided into five parts; 

1. Reliability analysis basic concepts. 

2. Formulation of the problem of determining system 
behavior by correlation of individual system component 
behavior. 

3. Mathematical determination of system reliability. 

4. Effect of preventive maintenance. 

5. The problem of interacting components. 
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PROCEDURE 

I. Reliability Analysis Basic Concepts 
Reliability and Ms an Time to Falliire 

Before introducing the method by which to solve the 
reliability analysis of a complex propulsion system, a 
short explanation of the theory behind the mathematical 
derivation of reliability will be given. This Is deemed 
necessary since reliability analysis is essentially a new 
field and although it has as a foundation the mathematics 
of probability and of statistics, there are some intrinsic 
differences. The standard reliability terms are defined in 
the Nomenclature section of this thesis, however since the 
basic mathematical formulation is short, it will be included 
in the text. 

To determine reliability mathematically, the following 
is presentedj 
Let 

X = f (x^,X2^ . . . be a vector of performance 

characteristics of the 
component 

A = region of the sample space representing the 
satisfactory performance of the component, 
w = f (w^,W2, . . .,Wj^)\ be vectors specifying the 

> upper and lower limits on 

z = f (z^,Z2^ o . o,Zj^) J gatisfactory performance 
t = time of the interval under investigation 

f(x,t) = probability density function of the performance 

characteristics 

R = reliabilltv of the comnonent 
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therefore 

R(x C Aj t) = Prob.(w < X < z) = / f(x,t) dx 

” A 

employing the proper definitions [ 2 . 0 ] 

R("t) = f(x) dx = 1 - /^ f(x) dx = 1 - F(t) 
t 0 

furthermore _jt ^ 

R(t) = e ° 

where 

R^tl “ “ dF R(t))o 

-/^ p(x) dx 
f(t) = p(t)e ® 

The other parameter of reliability mathematics that is 
of great importance is the components' mean time before 
failure (MTBP) . 

If T^= mean time before failure 

E = expected value 

then 

= E(t) = /” t f(t) dt = t d^(t) = f°° R(t) dt 
^ 0 0 0 

In the event that the component has been operating for a 

period previous to the time of investigation; 

Defining t = 0 time component is new^ or system starts 

operating 

t = t^ time of investigation 
t = time to point of interest 
R(t/t^) = reliability of the component at time t given 
the component has operated to time t^. 



thus 
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and 




X R(t) 



f(tj dt 



= nfiT r " 



Component Failure rdstritutions 

Prom the preceding paragraphs it can he seen that relia- 
hillty and mean time to failure can be evaluated from the 
probability density life dlstrlbjition of the component. In 
the following paragraphs some of the more Important mathe- 
matical approximations of component life distributions are 
discussed. These distributions have been derived from ex- 
tensive analysis of dperating data and life tests. It must 
be remembered^ however^ that they are just statistical models 
characterizing a physical phenomenon and are, at least, a 
good approximation of the reality. For a better understanding 
of the real situation, confidence limits are associated with 
each of these approximations. 

Normally, the characteristic parameter used to describe 
life distributions is the instant failure rate, p(t). The 
time dependence of p(t) that is representative of a normal 
component is shown in Pig, I. 
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FIGURE I 

Representative p(t) vs t plot 




where 

I = early failure region which results from manufacturing 
flaws, In-translt damage, etc. which cause normal 
failure at an early stage In the operating cycle. 
Usually In real situations there Is a "debugging" or 
run-in period which eliminates this problem. 

II = chance failure region caused by random failures which 
are normally considered Independent of time and 
having a fixed probability of occurlng at any time In 
this region. 

Ill = wear - out region due to the gradual deterioration of 
the component with time. 

Exponential Distribution 

This Is the most widely used life distribution In 
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rellability analyses. The reasons for this are It," 

1. characterizes components that have been optimized 
to limit failiires, 

2. characterizes components having a predominance of 
hmnan errors, 

3. characterizes complex components, 
characterizes components consisting of parts of 
mixed ages, 

5. approximates what is usually the behavior in the 
chance failure region, 

6. is the easiest mathematically to work with. 

The equations associated with the exponential distribution 
arej 



p(t) = constant = c 
f(t) = p(t) 

R(t) = 

^f = pW 



FIGURE II 



Exponential Distribution 







This distribution is generally employed to describe the 
wear-out phenomenon. It also has been found to apply to; 






1. components exhibiting homogeneous deterioration 
properties, 

2. components subjected to small variations in environ- 
mental severity, 

3. components whose failures occur at times well 
removed from t = 0, and whose MTBF is large 
compared to its standard deviation. 

The equations associated with the normal (Gaussian) distri- 
bution are I 



M = mean lifetime 



a = standard deviation of M 



e 




p(t) 







t 






cr*s/^ t 



= M 



if erf (x) = § e“^ du 

» r\ 



0 



then 



R(t) = I [1 - erf(£^)] 
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PIGURE III 



Normal Distribution 






M t t t 

Welbull Distribution 

Although this distribution has rarely been used in 
reliability analyses to date, statistical investigation of 
life tests of mechanical tests shows that this distribution 
seems to be applicable in many cases. The distribution can 
be described by a two parameter method or a three parameter 
method. It includes as a special case the exponential 
distribution and can approximate a normal distribution 





O = 3.25). 

For the three parameter casei 

where a = scale parameter ^ ^ 

p = slope parameter 
y = location parameter 



p(t) = 
f(t) = 



e 




5 



- 20 - 



it-yf 

''a 

R(t) = e 

= 7 + r(p+l) 



FIGURE IV 

Welbull Distribution 




Some other distributions occasionally used in relia- 
bility studies are I 

1 . lognormal 

2. gamma or beta 

3 . mixed Welbull 

however, the exponential, normal, and Welbull are the most 
Important of the distributions and are the ones used for 
examples in this thesis. 



I 
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II Formulation of the Problem 
Defining the Problem 

One of the most critical parts of performing an analysis 
of any system is the complete definition of the problem to be 
studied. This definition basically must answer three 
questions s 

1. What are the objectives of the analysis? 

2. What components are to be considered comprising the 
system? 

3. What will be considered the failure events of the 
system and components? 

These must be clearly and exactly defined if clear and 
exact answers are to be expected. In some cases if Just 
figure of merit or trend results are desired, greater liber- 
ties may be taken in the definitions. Once these three 
questions are answered, the problem has been defined and the 
results of the Investigation are completely dependent upon 
these definitions. 

Objectives 

When defining the objective of the analysis, the question 
is, what is wanted and what should be expected from the 
analysis. In order to better illustrate how to accomplish an 
operations analysis of a system, a hypothetical marine pro- 
pulsion plant will be postulated and investigated. It must 
be realized by the reader that this example will be necessarily 
simplified since the object is not mmierical results but . 
Instead the method of the solution. 
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The objective of this analysis will be to find the 
reliability of a combined diesel and gas turbine marine 
propulsion plant. The result desired is how does the 
reliability of the overall plant vary with time in oper- 
ation and how should one go about improving this relia- 
bility. Naturally in this example the value of the 
numerical results- are useless since there is no statisti- 
cal basis for the numbers used. In an actual analysis, 
statistical life tests would have to be performed on all 
the components involved. From these statistical tests, 
confidence limits can be found and applied to the final 
results. Methods of finding confidence levels are dis- 
cussed in Lloyd and Lipovj [ ? ] . 

System 

The next step in the problem is to define the system. 
The system can be defined in as basic or complex terms as 
the investigator deems necessary in order to obtain the 
results desired. In this example the system will be com- 
prised of the following components ^ 

1. diesel engine 

2. lube oil pump 

3. fresh water pump 

4. salt water pump 

5 . heat exchanger 

6 . gas turbine 

7. reduction gear 

8. standby lube oil pump. 
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Prom the definition of reliability it is seen that in 
defining the system, the operating conditions which are 
within the prescribed limits must also be enumerated. 

There are five basic operating conditions that should be 
taken into accoiuitj 

1 o environmental 

2. conditions imposed by the operator or user of the 
system 

3. functional dependence of the components within the 
defined system 

4. interaction of component failure distributions 

5. maintenance policies. 

It is altogether possible that more than one combi- 
nation of these conditions may want to be investigated. In 
this case different reliabilities can be expected for the 
different combinations. For the purpose of this investi- 
gation the following operating conditions are defined; 

Environmental ; The only extraordinary environmental 
conditions imposed on the system will be the possi- 
bility of shock with the probability of of 

occurring. 

Imposed conditions ; The system must be able to 
operate for 20^ of its total operating cycle above 
cruising power. 

Functional dependence ; The functional dependence of 
system operation upon the components is; 

1. Diesel must perform successfully for the system to 
operate . 

P RAaiir-hlrvr» o*Aar> rmist ruerform succeasfullv for the 
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system to operate. 

3o Heat exchanger must perform successfully for the 
system to operate. 

4. Salt water pimip must nerform successfully for the 
heat exchanger to operate. 

5. Fresh water pump must perform successfully for the 
heat exchanger to operate. 

6. Gas turbine must perform successfully when more than 
cruising power is required. 

7. Lube oil has to be available to the diesel at all 
times . 

8. A defined shock level with a probability of P_ of 
occurring will cause failure of the reduction gear 
and heat exchanger. 

Interaction effects s No interaction among the component 
failure distributions will be assumed. This is not 
necessarily the case in the real system but at the 
present state of the art# interaction effects cannot be 
handled except by using "guess" factors. In a later 
part of the thesis, these Interaction effects will be 
investigated. 

Maintenance policies % For the propulsion plant of this 
example, no true maintenance actions will be considered 
performed. It will only be assumed that minor adjust- 
ments for power, temperature, and fuel regulation will 
be made. Of course in most systems, maintenance actions 
will improve the reliability and this effect will be 
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covered in another section of the thesis. 

Concept of Success and Failure 

In order to determine the failure distribution of a 
component, it is necessary to define what is considered a 
failure and what a success in the operation of the component. 
It is not necessarily true that the only failxire of a com- 
ponent is when it stops functioning because in many 
situations, the system will also fail if a component oper- 
ates outside certain tolerances. For instance, if a lube 
oil pump doesn't supply the oil at the required pressure 
even though the p\mip may still be running, the system will 
fail. It is therefore evident that the investigator must be 
careful in choosing the ground rules for success or failure 
of any test and relate them to the various modes of operation. 

To find the failure distributions relating to the com- 
ponents, a four step procedure must be followed; 

1. Define success and failure of the entire system ^ 
being analyzed for a particular mode of operation, 

2. Relate the definition determined from 1 to each 
individual components operation so that success -failure 
criteria may be ascertained for each component, 

3. Establish the type of statistical distribution which 
describes the failure phenomena. 

4. Estimate the parameters which completely define 
these distributions. 

After this is accomplished, the investigator is ready 
to begin the mathematical aspect of the analysis since all 
input data are now available. For the example of this thesis, 
it will be assumed that the above steps have been carried out 
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and the failure distributions of the components in question 
are knovni. The case where the out of tolerance distribution 
is independent of and different than the "catastrophic'' 
failure distribution will be taken into account in one case. 

Reliability Schematic Block Diagram 

The reliability schematic block diagram sometimes called 
the reliability structure model is a picture form of the 
functional relationship of the system and the components 
under investigation. It is not a necessary facet of an oper- 
ational analysis , however in complex situations the model 
will provide an insight into how to best attack the problem 
from the mathematical standpoint. Prom this diagram it is 
also easier to recognize the weaknesses or potential weak- 
nesses of the system from a reliability sense. The relia- 
bility schematic block diagram for the thesis ex^ple is 
shown in figure {Y, . 

FIGURE V 




I 
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III Mathematical Procedure 

Three Basic Rules of Reliability Mathematics 

Before the mathematical analysis of the system Is begun, 

a statement of the basic reliability rules of probability 

calculus will be discussed. 

The first rule states that If there are N mutually 

Independent components, each having a reliability of Rj^C'^), 

and that for satisfactory operation of the system all the 

components must function properly; then the overall system 

reliability R (t) Is 
sys 

N 

R (t) = TT Ri(t) 

sys 1=1 ^ 

This Is known as the Product Rule of Reliability and a 
system of this type Is called a serial system. 

The second rule states that If there are N mutually 
Independent components each having a reliability of Rj^(t) and 
unreliability of Qj^(t) so that 

Ri(t) + Q^(t) = 1 

and that for satisfactory operation of the system, only one 

of these components must function properly, then the overall 

system reliability R (t) can be determined since 

sys 

N JL. 

Q (t) = TT ei(t) = TT (i-Ei(t)) 

sys 1=1 1=1 

R (t) = 1 - Q (t) 
sys sys . 
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N 

R (t) = 1 - TT (l-Ri(t)) 

sys i=l ^ 

This is known as the Product Rule of Unreliabilities and 
a system of this type is called a parallel system. This is 
usually described by the term redundancy^ relating to the fact 
that there are alternate components to help the system 
operate successfully in case of failure of one or more of 
the other components. The normal method for increasing 
system reliability is to add redimdancy to the original system. 

T?he third and most important rule is based on Bayes ' 
probability theorem. It states that the probability of 
system failure is the probability of system failure given 
that the component is bad times the probability the component 
is bad plus the probability of system failure given that the 
component is good times the probability the component is good, 

Qs(t) = Qg(t)|R^(t) X R^(t) + Qg(t)|Q^(t) X Q^(t) 

This is called the Conditional Rule of Reliability and 
can be applied to any system Including serial and parallel. 

Mathematical Analysis of the Complex System 

With these rules the analysis of the complex system can 
be begim. They will not account for every situation en- 
coimtered in the investigation but they will be the foxmdation 
for all the computational work done. 

For the purpose of simplicity, the components will be 
referred to by the number in the upper left hand side of the 
block in Figure 7 , l.e. 1 = diesel, 2 = gas txirbine, etc. 

It is seen that blocks A, B, C, D form a serial system such 
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that If any one of these fall the system falls. If the 
reliability of each of these four blocks could be found then 

R (t) = R,(t).R^(t).R (t).Rj^(t) I 

sys ^ ^ 

This will be the approach to detem Inlng the overall system 
reliability. 

Reduction Gear 

The first component to be examined will be the reduction 
gear. It Is known that the reliability of the gear Is R]^(t) 
under normal operating conditions. However, It Is also known 
that some predefined shock amplitude will cause the gear to 
fall. 

Define as the probability that a shock of the 

defined amplitude will occur 

1 - fs 

Reliability of D = Reliability of 4 given no shock 

occurs times the probability of no 
shock plus the reliability of 4 
given the shock occurs times the 
probability that shock occurs. 

- (R4(t)|P3 = 0) Qg + (R4(t)| P,^)P^ 

(Ri,(t)|P3 = 0) = R,,(t) Ri,(t)|Pg = 0 

therefore 

Rp(t) = R4(t) Qg II 

Heat Exchanger 

In finding the total reliability of block C, we again 
use the Bayes Theory approach of Rule 3. 
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\ R 3 (t)(if 7 is good)fi^(t)+R 2 (t)(if 7 is bad)-Q^(t) 



which leaves R 2 (t)(if 7 is good) to be found 
R 2 (t)(if 7 is good) = R 2 (t)(if 8 is good)-Rg(t) + 

R 2 (t)(if 8 is bad)-Qg(t) 

Rg(t)(if 7 is good) = R 2 (t)(if 8 is good)^Rg(t) 

as R 2 (t)(if 8 is bad) = 0 

but 

Rg ( t ) ( if 7 and 8 are good ) = R(j ( t )f P.g = o] 

.% Rc(t)(Pg = 0)= R3(t)R^(t)-Rg(t) 

and it is seen that the dependence of the heat exchanger on 
the fresh water and salt water pumps is the same as a 
subsystem composed of these three components in a serial 
arrangement. However, the shock hasn't been taken into 
account therefore 



Gas Turbine 

In the analysis of the block containing the gas turbine 
the first unusual case of mathematical reliability prediction 
appears where the three basic rules are not sufficient 
enough to give an answer. In this case, probability theory 
must be used to formulate a new type of reliability equation. 
This is due to the fact that the gas turbine is only used 
periodically and only needs to operate properly when called 
into use. 



= R 2 (t)(lf 7 is good)'R^(t) 
as R 2 (t)(if 7 is bad) = 0 




III 



-31- 



The reliability of the gas turbine can therefore be 
found by the following procedure j 
Define 

fgCt) = failure density function of the gas turbine 
^ 2 ^^^ “ probability that the gas turbine is not in 
use at time ^ 

U 2 (t) = distribution function of periods of non-use of 
the gas turbine 

Q 2 (t) = unreliability of the gas turbine. 

The situation can be described by the statement; a 
failure of the gas turbine at time 1 t causes system 
failure by t under one of two conditions: 

1. The gas turbine is in use at time t when it falls. 

2. The gasturbine is not in use at t but is called 
into use after t and at or before t. 

Condition 1 is characterized by 

Q^(t) = /^[l-P2(T)]f2('r)dT 
o 

Condition 2 is 

Q2(t) = F2(T)U2(t-T)f2(T)dT 

and 

%{t) = Q^(t) + Q^(t) IV 

Since the periods of non use are considered to be 
randomly spaced. 

Defining 

N 2 = mean period of gas turbine non-use 
the non-use density fimctlon is 
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and 



■f" n ” 

UgCt) = /^ iT ® ^ 

o ^^2 

UgCt) = 1 - e" 




o 



KLesel Engine 

In deriving the reliability equation for block A, 
encompassing the diesel engine, another situation arises when 
the three basic reliability equations are not adequate. This 
situation is caused by the lube oil-standby lube oil pump 
combination and this combination is referred to as a standby 
system. To find the reliability of such a system 
Define 



f^(t) = failure density function of the lube oil pump 
fg(t) = failure density function of the emergency lube 
oil pvunp 

Since the standby lube oil pump only operates at some time 
t^, when the lube oil pump fails, the true failure density 
functions for the system are 
f^(t^) lube oil pump 
fg(t-ti) = standby lube oil piunp 

As the only time the system fails is when both pumps fail, 
the system unreliability is fooand by Rule 2, 

« (t) - <t(t)Qg(t) = /* f fg(T-t^)dT 

standby ^ o ^ o 

and the failure distribution function of the emergency lube 
oil pump 
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R(t) = 1 - Q(t) 
consequently 

R (t) = 1 - /^ F^(t-tJf,(tT)dt_ V 

standby o ° 1 5 1 1 

Using the same type analysis as for the heat exchanges 

R.(t)=R^(t)R (t) VI 

standby 

Overall System Reliability 

To find the overall system reliability, equation I is 
applied 

R (t) = R.(t)*R„(t).R.(t).R^(t) 1(a) 

where 

R.(t) = RT(t)*R (t) 

standby 

RgCt) = [1 - 

R^(t) = R3(t)-R^(t)*Rg(t).Qg 

Rp(t) = 

Equation 1(a) Is applicable for the overall system 
reliability for any type failure distributions &>v the com- 
ponents t.e. exponential, normal, Welbull, etc. In most 
Instances, however, the solution of such complex equations 
necessitates the use of a digital computer. For the purpose 
of easy understanding and readily obtainable results, the 
failure distributions of the components will be assmed expo- 
nential. This is the standard method In the majority of 
reliability analysis and for most complex systems not too 
unrealistic, at least for a first approximation. The 
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assmptlon of the exponential distribution requires two basic 
postulates; 

1, The components are not operated long enough so that 
they begin to wear out due to aging. Prom a reliability 
standpoint this is a logical way to operate the components 
as their reliability greatly decreases as they enter the 
wear out regime. 

2. Any failure of a component is completely random, 
i.e. it is not time dependent. 

Using the exponential distribution, the following 
equations are arrived at; 






-c.t -c_t c, 

= e fe » + U 

o 5 



- e ^ )] 



See Appendix [ B ] for proof . 



Rg(t) = e ^^^[1+ P 2 (t)( i )3 + e ^ [P 2 (t)(j — )] 






^ " ^2 



See Appendix [B ] for proof, 



-(Co+C7+Co)t 
Rc(t) = ' 3 7 8, ] 



-Cii t 

RpCt) = ’ ) 



thus 
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^ -(C-+C-+C|.+C„+Co)t -c^-t 

R (t) = Q?e 1 3 4 7 8 5 



sys 



times 



Cgt P2(t)C2 



[e ^ (1+ 



^ - C 
N 2 ^2 



)] + e 



' 6 "S 

t 

N 2 ^P2(t)C^. 'j 



"^6^. 



1_ 

N. 






Quantitative Analysis 

In order to more graphically show the value of the relia- 
bility analysis 5 representative values for the component 
instant failure rates (p2(t) = C^) will be chosen. The 
overall system reliability as a function of time can then be 
evaluated. 

The values chosen are; 

Cg 

^3 



'5 

'6 

7 

"8 



.0002 


hr”^ 


.0003 


hr -1 


.00006 


hr“^ 


.00004 


hr“^ 


.0004 


hr“^ 


.0005 


hr“^ 


.0005 


hr -1 


.0005 


hr“^ 


Co 

00 




0 

00 

• 




200 hr 




assume 


the 



into accoimt. This component's failure due to wear-out shall 
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be described by a normal distribution with a mean life time 
of 1200 hrs and a variance of 400 hrs. 

Figure VI shows R (t) vs t for the cases with and with- 



out one component wear-out consideration. 

For the purpose of understanding the individual com- 
ponents relationship to the system as a whole, each of four 
component failure distributions have been varied to find its 
total effect on the system. These four components were chosen 
due to the different nature of their operation; 



sys 



Component 

Diesel 



Time Operation 

continuous, dependent upon L. 0. system 



Gas turbine periodic 



L. 0. pump continuous with standby system 
F. W. pump continuous, nondependent 
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IV Preventive Maintenance ; 

In most of the reliability work done In the past, little 
Interest has been given to the effect of preventive mainte- 
nance . The reason for this Is twofold due to the fact that 
most of the reliability study has been applied to the missile 
and electronic fields 

1. In the missile field the operation Is usually "one 
shot" l.e. the component or system having once begun Its 
operation cannot be maintained, and after terminating this 
operation It Is unusable. 

2. In the electronic field, the failure distributions 
of the components are considered exponential thus making 
preventive maintenance useless. This Is so because the 
Instant failure rate Is Independent of time and therefore 
replacing the component at some specific time when It Is 
still operating successfully has no advantageous effect from 
a reliability standpoint. 

These two conditions do not necessarily hold however, 
for a marine propulsion plant. Such plants definitely are 
not "one-mlsslon" devices and It Is doubtful If the exponen- 
tial distribution will be applicable to the mechanical com- 
ponents found In such systems. It Is proven by R. F, Drenlck 
[•^ ] that If the Instant failure rate Is constant or decreases 
with time, then preventive maintenance Is not beneficial. In 
mechanical components. It. has been determined that the Instant 
falliire rate normally Increases with time. In this case pre- 
ventive maintenance can be worthwhile and should be considered 



! 
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in any operational analysis of a marine propulsion system, 

A maintenance investigation should provide two basic 
parameters i how valuable to system performance is preventive 
maintenance, and at what time intervals should it be per- 
formed. For the purpose of analysis two different replacement 
policies will be explored. The mathematical analysis is from 
the work of R. Barlow and L, Hunter [ 5" ] . 

Policy li 

After every t^ hours of continuing operation without 
failure, preventive maintenance is performed. If there is a 
failure before t^, the maintenance is performed at this 
time, t^, and the next maintenance is scheduled at time 

*1 + *o- 

Defining 

= expected fractional amount of time the system 

is operating as t — > » 

1 

the two parameters that determine the condition of the problem 
are and t^. The criterion of optimality will be the 

maxlmlgatlon of . 

T = expected time needed to perform repair to the 
8 

system after a failure 

T = expected time needed to perform a scheduled 
preventive maintenance action. 

By making the asstunption that after each repair action, the 
system is as good as new, the following equations are appli- 
cable for optimal results ^ 
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p(t ) A R(t)dt - P(t ) = 

o e s 

> T„ 

1 © s 

’’l” " l+(Te-Tg)plto) 

if T < T . no preventive maintenance should be made and 

© o 



Policy II ; 

In this case, preventive maintenance is performed on the 
system after it has been operating for a period of t° hours 
regardless of the number of intervening failures. If there 
is a failure before t° hours have elapsed only minimal 
repair is made which in effect puts the system back in oper- 
ating condition but does not change the basic system failure 
rate. 

Using the same definitions for t) and T and defining 
^m “ expected time needed to perform minimal repair 
to the system 

the following equations are determined 

t® T 

/ tp'(t)dt= 
o m 

P(*°) 

The question naturally arises as to which policy should 
be used for a particular situation. This question can be 
resolved by equating = r\ 2 ^ for a given T^ which is 

known and the same for both policies. By substitution of 



<oo 



w 
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values for T , T can be determined. By graphing 
Tg vs for a given which policy to use becomes 

evident. Figure XI shows such a graph for a Weibull 
distribution. 

A Maintenance Policy by Dynamic Programming 

Assume that a preventive maintenance policy is to be 
determined for a system having an instant failure rate that 
is a linear fimctlon of its operating time. The system is 
to be operated in discrete time intervals of length t^. 

If the system is repaired between operating periods the cost 
associated with the repair is C^. If the system fails 
during an operating interval it cost C^. for the repair and 
the system cannot be operated again till the next interval. 
Define 

C = cost of scheduled maintenance 
r 

= cost of "in service" repair 
^m ~ Probability of failure in Interval m + 1 
given it has not failed in m Intervals 
since the last maintenance 
f^(m) = least expected cost of making n more 
intervals given m Intervals have been 
made. 

There are two policies that can be followed, either 
repair the system at m or do not repair it. By the 
methods of dynamic programming covered by R. Bellman [ 6 ] 
the following dynamic programming equation is developed^ 
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f (m) = min 
n' 



'R; repair 



__NR: non-repair 

C^ > Cj, for logical results otherwise preventive mainte- 
nance is not worthwhile. 

To show mathematically the effect of the failure repair 
cost, C«, and the scheduled maintenance cost, C_, on the 
maintenance policy to be followed, the following problem will 
be postulated. 

A system has a failiire distribution curve that can be 
described by the Weibull approximation with the following 
parameters," 



y = location parameter = 0 
P = shape parameter = 2 

4 

a = scale parameter = 57 x 10 hrs 

= characteristic lifetime = r\ = 755 hr 
M(t) = N(t) = interval time = 100 hr 



therefore 



P(t) = 

57(10^) 



R(t) = e 57(10“) 

The question to be answered is how does the maintenance 
policy vary as a function of m, n, C|., and C^, Figures 
XII through XVI give the results of this investigation. 
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V Interaction of Component Failure Distributions 

In the previous sections of the thesis, the overall 
reliability of a hypothetical complex marine propulsion 
plant has been derived by the proper mathematical combi- 
nation of Individual component reliabilities. In actual 
situations however. It has been found that although this 
method Is mathematically correct. It does not yield the 
actual reliability observed In practice. Intuitively, It 
might at first glance appear that this poor correlation Is 
because the model Is not a good functional representation 
of the real system. Although this Is a possible reason, 
further analysis may show that mathematically all work Is 
accurate. Where, then. Is the mistake In the analysis? 

One of the major possible explanations and one that although 
recognized Is not Included In present reliability studies. 

Is that when components are Integrated Into the system, 
there Is an Interaction among them which can change their 
failure distribution. Because of this, the failure distri- 
bution observed In the single component life test Is not 
applicable and must be modified. The problem of finding what 
type of function this modification factor Is, will In actual 
practice be very difficult. However, If reliability Investi- 
gations are to be accurate, component failure rate Inter- 
action has to be taken Into account as It may have a 
pronounced effect on system reliability, preventive mainte- 
nance policies, and system availability calculations. It Is 
therefore Imperative that laboratory life tests of components 
suspected of Interaction and datft assimilation of real system 
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operatlon be analyzed so that interaction fxmctlons can be 
determined . 

For the purpose of illustration, two possible situations 
where interaction may play a substantial role in an investi- 
gation are given. 

a. The failure distribution of a diesel cylinder liner 
will be affected by the wear of the jsiston rings, 
the wear of the connecting rod bearings, and the 
fuel quality. 

b. The failure distribution of a reduction gear will 

be dependent upon the wear of the shaft bearings and 
coupling operation. 

Since no data is available on interaction, for the 
purpose of investigation only hypothetical Interaction 
functions can be postulated. In the following sections such 
assimied fvinctlons are examined. 

Effect of Interaction on Component Reliability 

The time dependent reliability of a component is de- 
scribed by 

R(t) = exp[- p(t) dT] 
o 

which is a function of a single variable, p(t) = the Instant 
failure rate. An equation Illustrating the main factors that 
influence component failure rate is characterized by 
p(t) = f(a, p, (|), t) 
a = derating factor 
P = environmental stress factor 
|) = interaction factor 
t = operating time 
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where 

a = accotints for the fact that the component can operate 
at different outputs 

P = is a function of the operational environment the 
component experiences, i.e. temperature, vibration, 
and acceleration stress levels 
(() = accovints for the effect of interaction among the 
component failure distributions 
t = time the component is In operation 

It is normally possible to find the Instant failure rate 
as a function of derating factor, stress factor, and time 
from experimental life testing of the component. An example 
of such a function can be shown in a diagram for a component 
that follows the exponential failure distribution. 




j 
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Because reliability engineering is relatively new, little 
experimental data is available to find such functional re- 
lationships and what is available only applies to the expo- 
nential failure distribution* It is valid to assume that the 
same type of relationships or equations can be foimd for other 
types of components and failure dsltrlbutions. The problem 
still presents itself •' how can the interaction factor be 
handled? Since no work has been done on this matter, only 
conjectural relationships can be assioned for the purposes of 
investigation. Such relationships will therefore be assiimed 
and their properties analyzed. 

Define 

A(t) = component Instant failure rate without taking 
into account this interaction factor, (|). 

therefore 

p(t) = f ( A(t), (|)) 

It is this function that will be hypothesized. In all 
eases X(t) will be considered a constant for ease of mathe- 
matical computation. The effect of the interaction factor 
which is determined in this exponential case will apply in 
effect to other failure distributions. 

Investigation 1^: <|) = additive and constant 

p(t) = A + 2 C. 

n ^ 

R(t) = exp[-(^ + 2 C^)t] 

n 

T^ = 1/(a + 2 C. ) 

n 
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Investlgatlon 2: (j) = additive and time dependent 

p(t) = A + 2 C. t 
n 

R(t) = exp[-^t] exp[- i 2 C.t^] 

n ^ 



- exp[g 2 z j 

See Appendix C bJ for proof. ” 

Investigation 3.: <}) = multiplicative and time dependent 



P(t) 

R(t) 



Tf = 



= exp[- 
1.254 






Investigation (|) = additive and a function of a constant 

and a time dependent term 
p(t) = A + + Cgt 

from correlation with investigation 2 

R(t) = exp[- (A+C^)t] exp[- ^ Cgt^] 

T, = expti^l (1 - 

It can be seen that these values for R(t) and 
can differ greatly from the values obtained for the standard 
exponential failure distribution, l.e. 

R(t) = exp[-?^t] 



Figiire XVII shows the relative importance of 2 C. for 

n 

investigation 2 where A = .01. 
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Preventlve Maintenance Policy by 
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figure XVI 

Preventive Maintenance Policy by 
Dynamic Programming for ratio 
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DISCUSSION OF RESULTS 
Reliability Analysis 

In Figure VI the overall system reliability as a function 
of time Is shown. The plot graphically displays what degree 
of reliability can be expected from the system assuming; 

a. the values for the Instant failure rates specified 
previously. 

b. the components are good as new at t = 0. 

c. no preventive maintenance Is performed on the system 
during the operating cycle. 

It can be determined that this overall reliability curve 
for the system consisting of different components with unlike 
operating conditions can be approximated by; 

R (t) = e“*^^ 
sys 

where C = .OOI98 

Other results that can be obtained from Figure VI are; 

1. If the Instant failure rates of all the components 
were decreased by 50^, the overall system reliability will be 
Increased greatly as ^ Increases, l.e. 

t = 200; Rgyg is Increased by 20.2^ 

t = 50O; I^gyg is increased by 61.2^ 

t = 800; Rgyg is increased by ll6.8^ 

2. If one component has wear-out characteristics (fresh 
water pump) such that It should be Included In the reliability 
analysis, the overall system reliability Is adversely affected. 
This effect becomes especially prominent as the system 
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operating time approaches the mean life time of the component 
and Is also dependent upon the varlence of the mean life time. 
Figure VI shows that a"t t = 800, Hgyg is decreased by 15.7^ 
due to wear-out. The equation that Is applicable when Investi- 
gating a system or component which exhibits random and wear- 
out failures during the operating cycle Is; 

R.(t) = R. (t)-R. (t) 
w r 

defining Rj^(t) = overall reliability of component 1 

(t) = reliability of component 1 exhibiting 
^w 

only wear -out failures 

R. (t) = reliability of component 1 exhibiting 
^r 

only random failures 

Since the wear-out reliability factor Is multiplicative, 
the system reliability will decrease greatly with time as the 
number of components experiencing wear-out Increases. 

3. Figures VII through X Illustrate the relative 
Importance of each component to overall system operation 
from a reliability viewpoint. These figures show percentage- 
wise how the decrease In a specific components Instant failure 
rate will affect the system by using the formula; 



R = sys 
sys 



^sys^"^^’^ decreased p(t))-Rgyg(pi(t)=C^) 



R3y3(Pl(t)=Ui) 



It Is observed that a decrease In the fresh water pumps 
failure rate by 50^ Increases the overall system reliability 
by 22^ for an operating cycle of 800 hours while the same 
decrease for the diesel engines failure rate only Increases 
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the reliability by 7.8^. Therefore if the system reliability 
needs to be Improved, strictly from a reliability standpoint 
of the four components investigated the fresh water pump per- 
formance should be Improved first and then the lube oil pump, 
gas turbine, and diesel respectively. An investigation like 
this immediately indicates where in the system ^improvements 
should be made to increase overall reliability. Conversely, 
it also reveals what increase can be expected in system relia- 
bility for an improvement in a specific components Instant 
failure rate. Consequently, graphs of the type shown in 
Figures VI through X, which can be determined for any system 
if the instant failure rates of its components are known, 
provide an intimate insight into system behavior and provide 
valuable information for decisions regarding engineering 
Improvements . 

Preventive Maintenance 

Figure XI shows the relative merit of using either of 
the two policies discussed on pages -38 of this 

thesis. This graph is representative of the preventive 
maintenance situation for a system having a failure rate 
that increases with time. As an example of how to use the 
figure, asstunlng T_ = 7 hours and the expected value of 

o 

Tg = 4 hours then policy I is applicable imless Tj^ £ 2 
hours. Similarly, for any value of T there is a corre- 
spending value of T^^ which determines the boundary for 
choosing the optimum policy to follow. It will generally 
be found that policy II is most applicable to complex systems 
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whlch have many subsystems liable to fall and where keeping 
records of times to failure for each sjibsystem becomes 
cumbersome. It Is Interesting to note that If costs of 
falliire are substituted for the times of repair, then these 
same equations minimize the expected cost of operating the 
system. 

The method of determining preventive malfitenance 
scheduling by dynamic programming Is more suitable to a 
marine propulsion plant than the previously discussed method. 
The reason for this Is that with operating schedules to main- 
tain It Is probable that optimum preventive maintenance times 
fo\md will not coincide with the ships schedule. This method 
also provides a continuous picture of what policy to take re- 
garding preventive maintenance over all combinations of tlifie 
Intervals (trips) completed and time Intervals still to be 
completed. Figures XII through XVI Illustrate the results of 
such an analysis made for a system having the Welbull failure 
distribution described In the Procedure section. Prom these 
graphs the expected cost of maintenance can be minimized for 
any operational situation of this system If the cost of 
failure and of maintenance are known. These figures are only 
applicable for m, n > 0 as If m = 0, the system has Just 
been repaired and It does no good to repair It again, and If 
n = 0, there are no more trips to make so the expected cost 
= 0 . 

Although the graphs are self-explanatory, one Interesting 
observation can be made which Is not Intuitively obvious nor 
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would normally be expected if the investigation had not been 
made. Take as an example. Figure XIII for = 10 and the 
cxirve representing = 2, It is seen that if 8 trips have 
been made (m = 8) and 2 more are planned (n = 2), the policy 
should be to repair the system. However, with the same past 
history and if 3 more trips are necessary (n = 3)^ the system 
should not be repaired. The explanation of this Is; 

1, at large m the probability of failure in n inter- 
vals is relatively high and for small n, repair is the best 
policy since in effect it starts the system off at m = 0, 

2, as n Increases, this effect of starting at m = 0 
becomes less Important and the cost of preventive maintenance 
is the dominant factor. The reason is that at large n, even 
if maintenance is performed, there is still a significant 
probability of fallxire in the following n-1 intervals which 
causes a greater expected cost than if no repair is made and 

a high probability of system failure is accepted in the next 
n intervals , 

Analysis of Figures XII through XV provides the results 
for the determination of Figure XVI, Figure XVI shows that 
policy determination is solely dependent upon failure distri- 
bution, length of Inteirval, and the ratio of to C^. 

This fact lessens greatly the number of calculations that 
have to be made in a system analysis since assimiing the 
system failure distribution and Interval time are known, only 
one graph has to be derived instead of one for each possible 

value of C-, As the values of and are bound to 

I I r 




II 
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vary with time due to the normal trend of increasing repair 
and material costs, the advantage of the single graph repre- 
senting the ratio is evident. 

Interacting Components 

Although little if any investigation has been made into 
the problem of the interaction component failure distri- 
butions, Figure XVII Illustrates that in certain cases this 
effect can be very important. For instance, it is seen that 
if p(t) = .01 + 4(l0**^)t the mean time to failure is 

Tf 

decreased 20^. On the other hand for t = ^ , this inter- 
action effect only decreases reliability by less than 
Consequently although the effects of interaction may not 
affect the basic calculations for reliability significantly 
and therefore may be difficult to observe and determine, they 

f- 

definitely can become significant in the overall system 
analysis through which is used in many operations analy- 

sis calculations not discussed in this thesis, i.e. renewal 
theory. It therefore becomes imperative that investigations 
should be made determining any interaction present in a 
system if worthwhile results are to be obtained from an 
analysis of the system. 
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CONCLUSIONS 

In addition to the specific conclusions stated in the 
Discussion of Results, the following general conclusions can 
be drawn. 

1. The operational analysis of a marine propulsion 
plant will provide results not otherwise available which 
will be beneficial In making decisions regarding; 

a . preliminary design 

b. operational policies 

c . maintenance actions 

d. system Improvement and redesign 

2. In such an analysis the definition and formulation 
of the objectives of the Investigation determine all suc- 
ceeding stages of the analysis and the utility of the results, 

3. Statistical data are a necessary prerequisite to any 
operational analysis. The definition of the problem governs 
the areas of investigation for the statistical tests. 

4. Reliability as a function of time can be character- 
ized by the single parameter p(t), instant failure rate. 

The effect of an increase or decrease in this parameter is 
amplified as time increases, 

5. The wear-out failure distributions of the individu- 
al system components should be included in the analysis of 

a mechanical system as they can greatly affect the overall 
system reliability. This consequence becomes especially 
prominent as t approaches the MTBF of a component. 

6. Prom an operational analysis of a system, the 




j 
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relatlonship of each component's performance to overall 
system operation can be determined. Such information will 
be of considerable significance in the determination of 
policies regarding system design and improvement. 

7* Preventive maintenance is only advantageous for 
devices having an instant failure rate that increases with 
time, 

8. For this type device there is an optimal time for 
the performance of preventive maintenance which will maxi- 
mize availability. An optimal time may also be found which 
will minimize expected cost of operation, 

9. Dynamic programming can be utilized in determining 
an optimum maintenance policy for systems operating in 
discrete time intervals. This policy will give a complete 
picture at any time in the systems operating cycle of what 
maintenance action to follow in order to minimize the 
expected cost of repair, 

10. The effect of interacting component failure distri- 
butions normally neglected in operational analyses and 
difficult to determine quantitatively can considerably 
affect the results of the investigation and should therefore 



be included 




p 
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RECOMMENDATIONS 

1. The concept of the operational analysis of marine 
pro pills Ion plants should be accepted by naval engineers as 
an important tool in the evaluation of existing and proposed 
marine systems. 

2. A program of component testing and data assimilation 
should be initiated immediately for the purpose of determi- 
ning representative failure distributions of mechanical 
devices . 

3. Accepting the fact that corroborating statistical 
data are minimal, systems analyses should be carried out on 
existing systems for the object of developing and refining 
techniques and gaining facility in the procedures of such an 
investigation. 

4. Further examination should be made into the question 
of how out of tolerence l.e, wear-out failures are to be 
defined and handled. 

5. A detailed investigation should be made into the 
problem of failure distribution dependence. This should be 
explored both from a theoretically mathematical and an 
experimental viewpoint. 

6. The applicability of the method of dynamic program- 
ming to system analyses questions should be investigated. 

T 

Possible applications other than for maintenance policies 
might be; 

a, optimal reliability apportionment among the 
components , 
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b. optimal arrangement of components imder con- 
straints for maximm system reliability. 

c, optimal choice of components to minimize the 
deleterious effect of failure distribution 
interaction. 
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APPENDIX A 
Definitions 
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Fallure ; a detected cessation of ability to perform a speci- 
fied function or functions within previously established limits 
Independent Failure: those component failures which occur or 

can occur without being related to the malfunctioning of assocl 
ated components. 

Redimdancy : the existence of more than one means for ac- 

complishing a given task where all means must fail before 
there is an over-all failure to the system. 

Availability : the fraction the total desired operating time 

that a system is actually operable. 

Mean Time Before Failure (MTBF) : the mean or average time 

between successive failures of a component 

Cumulative Probability of Failure (Q(t), F(t)): probability 

that the life of a component is less than t. 

Probability Density Function of Failure (f (t) ) : probability 

that the life of a component will be between t and t + dt. 
Instant Failure Rate (p(t)): probability that a component 

will fail during the time Interval t to t + dt conditional 
upon it surviving up to time t. 

Reliability R(t) : probability that a component will operate 

within specified limits for the time, t, and operating con- 
ditions Imposed. 
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appendix B 
Details of Procedure 

I Periodically used component 

II Standby system 

III Interaction 

IV Computer program 
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I, Determination of R(t) for a component used periodically- 
assuming an exponential failure distribution. 



F(t) = [l-P(T)]f (x)dT + P(T)[G(t-x)]f (x)dT 

o o 



P(x) = P 

G(t-x) = 1 - e ^ 
f(x) = C e"^'^ 

Q^Ct) = [1-P]C e'^'^dx = [l-P][l-e“^^] 

o 



= 1 - P - e-« + P e-« 






P[l-e ]C e“^^dx 



= P - P e"^^[l + — ) 



IT 



t-C 



P(t) = 1 - e"^^[l + PCj-^)] + PCj-^) 



N 



IT 



- c 



IT 



- c 



R(t) = 1 - F(t) 



R(t) = e‘^^[l + )] - ) 



IT 



t-C t-C 



II. Determination of R(t) for a standby system assuming 
exponential failure distributions. 

p(t) = for main component 
p(t) = Cg for standby component 
For a standby system Q(t) = F 2 (t-t^)f^(t^)dt^ 
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. . -c^Ct-tJ 

PgCt-ti) = r f‘(t-t^)dt = r Cge dt 

tl 

-CT(t-t^) 

= 1 - e ^ 

Q(t) = r F2(t-t^)f^(t)dt^ = r [1-e ^ ]C^e ^ ^ 



dt. 



. -C^t^ -C^t . (Cp-CT)tT 

= r C^e ^ ^dt^ - C^e 2 r e ^ ^ ^dt^ 

o o 

t 



-C,t 



, , 1” (VCl)tl 

= 1 - e - M _A - \ e He 



] 



-C,t Ct -C,t -Cot 
= 1 - e ^ - yr-4r-[e ^ - e ^ ] 






R(t) = 1 - Q(t) 



-C^t Ct -Cgt 

e 



R(t) = e ^ [1 + TT-^] - 

Cg-Oi 



III, Determination of R(t) and for p(t) 

p(t) = X + C^t 

R(t) = exp[-/^(X + Cj^T)dT] 



X + C^t 



= exp[-/^ XdT - C.T dr] 
o o 

= exp[-Xt - T dT] 

o 

Q 

= exp[— Xt - t^] 

1 P ^2 

^-xt ? V 

= e e 



= /" R(t)dt 
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= exp[-(^t + i C.t^)]dt 
o 

From C.R.C. Standard Mathematical Tables [ 7 ] 

r exp[-aV]<te = = a8|L 

0 

it can be seen that 

1 



I + Xt = I C,[t 2 + t(+ ^ ® - ^ ")] 



= |c,[(t + ^)" - (^n 



■ I ■ ? 



therefore 



R(t) = exp[^^^] exp[- ^ Cj^(t + ] 

therefore 

T^ = /“ R(t)dt = exp[^^^] /" exp[- ^ C^(t + ]dt 



let 



X = t + ^ 

^i 

dx = dt 



t = 0 
t = 00 






X = 00 



/“ exp[- ■! C^(t +^) ]dt = /“ exp[- ^ C.(x^)]dx 
o ^i X 

which equals (dropping the subscript) 

X 

/” exp[^ x^]dx - exp[|- x^]dx 
o o 

from C.R.C. tables 



/“ exp[§ x^]dx = ^/|^ 
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now let t = ^/^ X 


X = 0 


o 

11 


dt = dx 


^ = C 


t = 


therefore 


A 




A 

f' exp[^ x^]dx = 
o 


*/ir 0 


exp[-t^]dt 


and 

~ exp[-t^ 

*/V 0 


]dt = erf[-^] 





thus 

re*p[fx2)dx=4j[l -erf(^)] 





(1 - erf[ 




and 




5“ > 









I 







w. 



I 



# 
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IV Main Body of Computer Program for Solution of Equation 

(fortran) 

DO 114 N = 1, N1 
JJ = 0 

DO ll4 M = 1, Ml 
IP(N-l) 100, 100, 101 

100 A1 = 0. 

GO TO 102 

101 J = N - 1 

A1 = C0ST1(J,1) 

102 A2 = CASREP + A1 
A3 = RATE(M) * A2 
IF(N-l) 103, 103, 104 

103 A4 = 0. 

GO TO 105 

104 K = M + 1 

A4 = C0ST1(J,K) 

105 A5 = l.-RATE(M) 

A6 = A5 * A4 
C0ST1(N,M) = A3 + A6 
IF(N-l) 106, 106 , 107 

106 B1 = 0. 

GO TO 108 

107 B1 = C0ST2(J,1) 

108 B2 = SKDREP + RATE(l) * (CASREP + Bl) 

IF(N-l) 109, 109, 110 
B3 = 0. 



109 
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go TO 111 

no B3 = C0ST2(J,2) 

111 B4 - l.-RATE(l) 

B5 = B4 * B3 
C0ST2(N,M) = B5 + B2 
IF(JJ) 112, 112, 113 

112 JJ = 1 
PRINT 20 
PRINT 21, N 

113 L = M - 1 

114 PRINT 22, L, C0ST1(N,M), L, C0ST2(N,M) 
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